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Abstract 

Existence  and  uniqueness  theorems  are  proved  for 
two  boundary  value  problems  for  the  axisymmetric  deforma- 
tion of  a  circular  membrane  subjected  to  normal  pressure. 
The  nonlinear  Ftippl  membrane  theory  is  employed.   The 
shooting  method  is  used  to  establish  these  results. 
It  is  also  shown  that  if  the  edge  is  free  to  move  in  the 
plane  of  the  membrane  then  a  necessary  and  sufficient 
condition  for  the  existence  of  a  unique  solution  is  that 
the  pressure  is  self-equilibrating.   The  shooting  method 
is  applied  to  obtain  a  numerical  solution  of  the  fixed 
edge  membrane  with  uniform  normal  pressure. 
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1.   Introduction. 

We  consider  two  boundary  value  problems  for  the 
axisymmetric  deformations  of  circular  membranes  due  to 
an  axisymmetric  pressure  applied  normal  to  the  plane  of 
the  membrane.   We  employ  the  Fiippl  [1]  nonlinear  membrane 
theory.   The  first  problem^  which  we  call  problem  H,  is 
for  the  fixed  edge.   This  is  the  problem  considered  by 
Hencky  [2]  for  uniform  pressure.   In  the  second  problem 
(Problem  S),  the  edge  is  restrained  from  deforming  normal 
to  the  plane  of  the  membrane  and  is  subject  to  a  uniform 
radial  tension,  S  >  0,  applied  normal  to  the  edge.   We 
prove  that  there  exists  a  unique  solution  of  each  of  these 
problems.   Furthermore  we  show  that  the  radial  stress  is  a 
monotone  decreasing  function.   If  S  =  0,  then  we  have 
Problem  0.   We  prove  that  Problem  0  has  a  unique  solution 
if  and  only  if  the  pressure  is  self-equilibrating. 

The  "shooting  method,"  is  used  to  establish  these 
results.   The  method  is  constructive  and  it  is  the  basis 
of  a  well  known  procedure  for  numerically  solving  two 
point  boundary  value  problems.   It  has  been  used  previously 
to  prove  a  corresponding  result  for  a  class  of  second  order 
boundary  value  problems  [3].   However  the  differential 
equation  for  the  circular  membrane  is  singular,  and  hence 
Keller's  theorem  [5],  does  not  apply.   In  addition  the 
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boundary  conditions  for  Problem  H  are  not  considered  in  [3]. 
Our  analysis,  which  differs  from  that  in  [^],    is  clearly 
applicable  to  a  class  of  singular,  two  point  boundary  value 
problems . 

Previous  existence  and  uniqueness  theorems  for 
Problems  H  and  S,  using  an  iterative  procedure  due  to 
Ghaplygin,  have  been  given  in  [4]  and  [5].   It  is  suggested 
in  [6]  that  corresponding  results  have  been  obtained  for 
Problem  0.   Dicky's  analysis  [7],  which  is  related  to  that 
in  a  previous  study  of  circular  plates   [8],  is  for 
Problems  H  and  S  with  uniform  pressure.   Furthermore,  for 
Problem  S  it  requires  that  S  >  Sq  >  0,  where  Sq  is  a  fixed 
niimber  depending  on  the  pressure.   A  simple  uniqueness  proof 
using  elementary  arguments  is  sketched  in  [9]. 

In  Section  2,    we  formulate  the  boundary  value 
problems.   Since  the  differential  equation  for  the  circular 
membrane  is  singular,  we  cannot  employ  the  standard  existence 
and  uniqueness  theorems  [10]  for  initial  value  problems. 
Thus  in  Section  5  we  establish  basic  theorems  for  the 
initial  value  problem.   The  main  results  for  the  boundary 
value  problems  then  follow  easily  by  the  shooting  method. 
They  are  proved  in  Section  4.   In  Section  5  we  apply  the 
shooting  method  to  obtain  a  niimerical  solution  of  Problem  H 
with  uniform  pressure. 
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2.   Formulation. 

A  circular  membrane  of  radius  R  and  thickness  t  is 
subjected  to  a  normal  pressure  p(r)j  where  r  is  the  radial 
coordinate.   We  assiime  that  the  deformation  is  axis ymme trie. 
This  implies  that  the  radial  and  circ\imferential  membrane 
stresses  a  and  a  and  the  radial  and  normal  displacements 

J.  Ty 

U  and  W  are  functions  of  r  only  and  that  the  circumferential 
displacement  and  the  shear  stress  vanish.   Then  Fttppl's 
membrane  equations  [1]  can  be  reduced  to  the  single  nonlinear 
ordinary  differential  equation  for  the  dimensionless  radial 
stress  s, 

(1)  Ls  =    (x^s'  )'  +  Q/s^  =  0  . 

Here  we  have  introduced  the  following  notation 

X  =  ^  ,      0<x<l^     s(x)  =    a^(r)/E 

Q(x)  ^  (Qo/x)[  J  (p(q/E)e  dif   >0,      Qq  ^  (2tV)-\ 

0 

where  E  is  Young's  modulus.   We  assxame  that  p(r)  is 

piecewise  continuous  on  [0,1].   Thus,  Q(x)  has  a  piecewise 

continuous  derivative  and  from  (2)  there  is  a  positive 

constant  H  such  that 

(3)  Q  <  Hx-^  ,    0  <  X  <  1  . 
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The  dlmensionless  normal  displacement  v  and  the  dlmensionless 
circumferential  stress  s   are  defined  by 

v(x)  H  W(r)/(2l/2j^),   s^(x)  =  aQ(r)/E  . 

They  are  determined  from  the  equilibriLun  equations 

(4)  V  ==  s-l(QA)l/2  ^     3^  ^  (^^),_ 

To  complete  the  formulation,  conditions  are  required 
at  the  center  x  =  0  and  the  edge  x  =  1  of  the  membrane. 
Symmetry  and  boundedness  of  the  stresses  and  displacements 
imply  that  s(0)  is  bounded  and  that 

(5)  s'(0)  =  0  . 

At  the  edge  we  consider  two  different  boundary  conditions! 

(6a)  s(l)  -  S  >  0 

(6b)  s'(l)  +  as(l)  =0,    a=l-V>|, 

where  v  is  Poisson's  ratio.   We  also  require  that  in  both 
cases  v(l)  =  0.   This  is  the  initial  condition  to  determine 
V  from  (4).   The  condition  (6b)  implies  that  the  radial 
displacement  vanishes  at  the  edge.   This  is  the  fixed  edge 
case  studied  by  Hencky  [2].   The  prescribed  constant  S  is 
proportional  to  the  radial  stress  applied  to  the  edge. 
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Thus  we  define  two  boundary  value  problems.   They 
consist  of  determining  functions  s(x)  with  continuous  first 
derivatives  on  the  closed  unit  interval  [0,1]  which  satisfy 
(1)  on  (0,1)  and  (5)  and,  for  Problem  S  condition  (6a), 
and  for  Problem  H  condition  (6b).   Since  s  satisfies  (l) 
on  (0,1)  and  Q  is  continuous,  s"  is  continuous  on  (0,1). 
The  requirement  that  s'(x)  is  continuous  on  [0,1]  insures 
the  boundedness  and  continuity  of  v  and  s   on  [0,1]. 

The  existence  of  a  unique  solution  of  these  boundary 
value  problems  will  now  be  demonstrated  by  using  the 
"shooting  method."   Thus  we  consider  the  Initial  value 
problem*.   Determine  a  function  w(x;c)  with  continuous 
second  derivatives  which  satisfies 

Lw  =  0  ,     x  >  0  , 
(7) 

w(o;c)  =c,    w'(o;c)=o. 

The  basic  idea  of  the  shooting  method  is  to  show  that  there 
exists  a  unique  value  c„  of  c  such  that 


(8)  w(1;Cq)  =   S  ,     for  Problem  S  , 


and  a  unique  value  c^  such  that 


(9)  w'(1;cq)  +  aw(i;cQ)  =  0  ,   for  Problem  H 
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Then  s(x)  =  vi{xlc^)   will  be  the  \inique  solution  of  the 
boundary  value  problem. 

To  apply  the  shooting  method  we  require  basic 
properties  of  the  solutions  of  the  initial  value  problem 
(7).   Since  the  differential  equation  is  singular,  we  cannot 
apply  standard  theorems  [10].   In  the  next  section  we  shall 
therefore  analyze  (7). 
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5.   The  Initial  Value  Problem. 

The  function  Q(x)  is  defined  in  (2.2)  for  x  in  the 
unit  interval.  It  is  convenient  to  extend  Q  continuously 
to  X  >^  1  by  the  definition 

R 

(1)  Q(x)  =  (Qq/x)[  J  ^pd^  +  x^-  1]2  ,   x>l. 

0 

Then  the  inequality  (2.3)  is  preserved  for  all  positive  x, 
i.e. 

(2)  Q(x)   <  Hx^  ,    X  _>  0  . 

It  can  be  shown  by  direct  calculation  that  (2.7) 
is  equivalent  to  the  integral  equation 

X 

(5a)  vj(x;c)  =  c  -  J  k(x,|)w"^  d?  , 

0 

(3b)    k(x,|)  =  I  Q(e)(r^-x-2)  >  0  ,   f or  0  <  e  <  X  , 

in  the  sense  that  every  solution  of  (2.7)  is  a  solution 
of  (3)  and  conversely.   By  differentiating  (3a)  we  obtain 

X 

(^)  W (x;c)   =   -  x"^  J  Qw"^  d^  . 

0 

Thus,  every  solution  of  (3)  is  a  monotone  decreasing 

function.   If  Q  =  0  in  [0,1],  then  from  (3),  w  =  c  is  the 

unique  solution  of  (2.7).   Then  Problem  S  has  the  unique 
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solution  s(x)  =  w(x;S)  =  S,    and  Problem  H  has  the 
unique  solution  s(x)  =  w(x;0)  =  0.   We  shall  hereafter 
assume  that  Q(x)  ^  0  in  [0,1]. 

The  basic  properties  of  the  initial  value  problem 
(2.7)  are  contained  in  the  Theorems  I-5. 

THEOREM  1.   For  every  S  >  0  and  every  c  >  S, 
there  is  a  unique  number  x„  >  0  such  that  only  one  solution 
of  (2.7)  exists  in  the  interval  0  j<  x  _<  x„  and  w(xq;c)  =  S. 
The  solution  is  a  positive  and  monotone  decreasing  function 
of  X  in  this  interval. 

ProofI   We  first  demonstrate  the  liniqueness  in  any 
interval  in  vjhich  w  does  not  vanish.   For  a  fixed  c,  let 
w,  and  Wp  be  two  distinct  non-vanishing  solutions  in  the 
interval  0  <_  x  <_  a,  where  a  >  0.   Since  c  >  0  it  is  always 
possible  to  determine  such  am  interval  by  making  a  suffi- 
ciently small.   Then  D(x)  =  w-,  -  w^  satisfies  the  linear 
integral  equation'. 


X 


(5)     D(x) 


k(x;^) 


0 


w 


l(?)+vv2(t) 


w^(l)   w2(^) 


D(^)  di   ,      0  ^  X  <  a. 


Since  the  kernel  is  continuous  in  the  triangle  0  _<  x  <^  a, 
0  j<  ^  _<  X,  the  Volterra  equation  (5)  has  a  unique  solution. 
Thus  D=0in0<x<a  and  the  initial  value  problem  has 


We  need  only  consider  c  >  0  since  for  c  j<  0  a  solution  of 
(2.7)  cannot  solve  the  boundary  value  problems. 
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no  more  than  one  solution  in  this  interval. 

■St 

To  prove  the  existence  we  consider  a  strictly- 


monotone  decreasing  sequence  of  initial  points  [a  }  with 
a-,  <  1  and   lim   {a  )  =  0.   For  each  n,    the  corresponding 


n  — >oo 
initial  value  problem  for  w 


n 


(6)    Lw^  =  0,  ^^i^n>^)   =  c,   %(an^°^  =  °  '      n=l,2,...  , 


^n-^-"n 


has,  by  standard  theorems  [10 ],  a  unique  solution  in  an 
interval  0  <  a   ^  ^  <_  a   in  which  w   does  not  vanish. 
Clearly,  (6)  and  the  integral  equation 

X 

(7)  w^(x;c)  =  c  -  /  k(x,|)w-2(e)  d?  , 

^n 
are  equivalent.   In  addition  we  have  by  differentiating  (7) 

X 

(8)  w^   =   -  x"^  J  Qw^^  dC  . 

a 
n 

If  Q(a  )  7/0  then  from  (8),  w   is  a  strictly  monotone 
decreasing  function  of  x.   Suppose  that  Q(x)  =   0  for  x  in 


An  existence  proof  can  also  be  obtained  by  showing  that  the 
initial  value  problem  has  a  convergent  power  series  expan- 
sion in  a  sufficiently  small  neighborhood  of  x  =  0.   Then 
choosing  an  initial  point  x  :=  x*  =/  0,    in  this  interval,  the 
coefficients  in  the  differential  equation  are  not  singular 
at  X*  and  standard  theorems  [10]  can  be  applied  for  x  >  x* . 
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the  interval  a  :!  x  ^  X  and  Q(x)  fi^  0  in  any  interval 


w 


ith  X  as  the  left  end  point.   Then  w^  =  c  in  a^  ^  x  _<  X^ 


and  it  is  strictly  monotone  decreasing  when  x  >  X  .   For 

every  positive  niunber  S  <  c  either  w  (x;c)  =  S  for  a 

unique  x  =  x   or  else   lim  w   =  0  with   lim  w„  =  A  _>  S. 

^         X  — >oo  X  — >oo 

The  latter  possibility  cannot  occur  since  by  (8) 

lim  w   =  -  ( — 5-)  lim  — ^  and  using  (1)   lim  w   =  -  oo  . 
x-xE  ^      3A^  x->oo  yi  X  ->oo 

Thus  in  (6)  and  (7),  a   >  x  . 

LEMMA  1.   The  sequence  [x  }  is  monotone  decreasing 
and  converges  to  Xq  >  0. 

Proof:   Suppose  the  sequence  is  not  monotone 

decreasing.   Then  for  some  n,  x   i  >  x  and  there  is  a 

point  X*  <  X  at  which  w  (x*  ;  c)  =  w  -^(x*;c).   This 
implies  using  (7)  that 


^n 


k(x%?)w;2^  di  =  J  k(xM)(w;;2-w;^^)  d^  <  0 


^n+l  ^n 

since  w^_^-^  j<  w^  in  a^  <_  x  _<  x*.   If  Q,  ^  0  in  (a^^^n+l^ 
then  this  is  a  contradiction  because  a   .  "^  ^n*   If  Q  -  0 

^^  (^n^^n+l)  ^^^^  ^n+1  "  %  ^^  ^^n^^n^  ^^^  ^n+1  =  ^n* 
Therefore  x   -,  <_  x  and  {x  3  is  monotone  decreasing. 

Since  x   >  0  for  all  n,  the  sequence  converges  to  x„  . 

The  limit  x^  must  be  positive.   This  follows  by  using  (7) 

in  w  (x  :c)  =  S  to  obtain 
n^  n-'  ' 
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X  X 

n  n 


(9)  0<c-S.  Jk(x^,^)w-2d^<-^J  ^^i^-^..-^2)  ai 

8S^   n   n 

We  have  used  w   >^  S  and  (2)  in  estimating  in  (9).   If  x„=  0 
then  the  right  side  can  be  made  arbitrarily  small  by  making 
n  sufficiently  large  since  a  — f  0.   This  is  impossible 
because  the  left  side  is  positive  and  the  lemma  is 
established. 

It  is  convenient  to  introduce  the  infinite  sequence 
of  functions  {t  (x))  defined  on  the  fundamental  interval 
0  j<  X  j<  Xq  for  n  =  1, 2,  ...  ^  by 


(10)  t^(x) 


c,       0  1  X  <  a^  ; 


w  (xlc),  a   <  X  <  x„   if  a^  <  x„ 
n^^'-'n—   —  0      n    0 


The  sequence  is  uniformly  bounded  by  c  since  w   is  monotone 
decreasing.   If  the  sequence  is  uniformly  Lipschitz 
continuous  with  respect  to  x  and  n  on  [0,Xq],    then  by 
Arzela's  theorem  there  exists  an  infinite  subsequence 
{t   (x)}  which  converges  uniformly  to  a  continuous  fiinction 
w„(x,c)  on  [OjX„].   The  function  w^  is  a  solution  w(x3c)  of 
(2.7).   To  see  this  let  n  ->  oo  over  the  subsequence.  Since 


The  sequence  is  uniformly  Lipschitz  continuous  with  respect 
to  X  and  n  if  the  Lipschitz  constant  K  is  independent  of  x 
and  n. 
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the  subsequence  converges  uniformly  and  the  integrand  in 
(7)  is  a  continuous  function  of  w^  .  we  can  interchange  the 
limit  and  integral.   Hence,  since  a^  — >  0,  w^  satisfies 
(3a)  and  is  a  solution  of  the  initial  value  problem*  (2.7) 
for  0  ^  X  _<  x„.   Furthermore;,  Wq{x'q)    =   S  and  the  mono- 
tonicity  is  a  consequence  of  (4). 

Therefore  to  prove  the  theorem  we  need  only  establish 
that  the  sequence  (lO)  is  uniformly  Lipschitz  continuous  in 
X  and  n.   This  is  easily  done  by  observing  from  (8)  that 
w   is  continuous  and  bounded  on  [a  ^x^].   For  any  x  in 
this  interval  we  have 

(11)  K\:-±,h'^^^±^i/^^^^i^"^^--' 

where  we  have  used  (2)  and  the  inequality  w^  ^  S,   Since 

w  =  0  for  0  <  X  <  a  ,  the  boimd  (11)  holds  for  all  x  in 
n  —   —  n^  ^   ' 

the  fundamental  interval.   Since  K  is  independent  of  x 
and  n  it  is  the  required  Lipschitz  constant. 

5f 

Every  convergent  subsequence  converges  to  a  solution  of 
(2.7)  in  0  ;<  X  ^  xq.   Since  (2.7)  has  a  unique  solution, 
the  entire  sequence  (lO)  converges  uniformly  to  the 
same  limit. 
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THEOREM  2.   Let  w(x;c)  be  the  solution  of  the 
initial  value  problem  (2.7).   Then  for  every  S  >  0  and 
every  c  >  S,  w(x;c)  and  w'(xjc)  are  continuous  functions 
of  c  for  every  x  in  the  fundamental  interval  [OjX„(c)]. 

Proof:   Consider  any  two  values  of  c  >  S,  c,  and  c 
and  let  Wj|^,(x)  =   w(x;c-j^)  and  ^^(x)  =   w(x;Cp).   Then  using 
(5),  the  difference  v(x)  =  w-,-Wp  satisfies  the  integral 
equation 

(12)  v(x)  =  c-^-c^  -  J  k(x;e)    -^2 2^ ^(^)  ^^  ' 

for  0  ;<  x  ^  m  =  min  [Xq(c-j^)  ,XQ(Cp)  ] .   Taking  absolute 
values  in  (12)  and  using  (2)  and  the  inequalities 
w^  <_  c^,    w^  ^  S,  i  =  1^2,   to  estimate  the  integrand 
we  obtain 

X 

(13)  |v|  <  Ic^-  Cgl  +  2A  J  ^\v\    ai,         2A  =  H(c-^  +C2)/(2S^) 

0 

Applying   Gronwall's    inequality    [10]    to    (I5)   yields 


(14)  Ivl      <      Ic^-cpl    e^^^ 


The  continuity  of  w  as  a  function  of  c  for  fixed  x  is  a 
consequence  of  (l4).   The  continuity  of  w'  follows  directly 
from  the  integral  relation  (4)  and  the  continuity  of  w. 
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The  final  three  theorems  are  concerned  with  the 

variation  of  the  point  x^(c).   We  show  that  x„(c)  is 

continuous  and  strictly  monotone  increasing  for  c  >_  S 

and  Xp,(oo)  =  oo  .   Clearly  if  Q(0)  4   0,   lim  Xp^(c)  =  0. 

If  Q(x)  =  0  for  X  in  the  interval  0  _<  x  ^  X  <  1  and 

Q,(x)  ^   0  for  any  interval  with  X  as  the  left  end  pointy 

then  lim  x  (c)  =  X  <  1. 
c^S  ^ 

THEOREM  5.  '^rX^)    is  strictly  monotone  increasing 

for  c  ^  S. 

Proof:   Assume  the  contrary^  i.e.  there  are  two 
values  C-,  and  Cp  of  c  for  which  c^  >  c-,  >_  S  and 
x„(c-,)  ^x:„(Cp).   Then  there  must  he  a  value  x  in  the 
interval  0  <  x  <_  x  (Cp)  such  that 

(15a)  w(x*;c-l)  =  w(x*;c2)  , 

(15b)  w(x;c-^)   <  w(x;c2)  ,  0<x<x*. 

Substituting  the  integral  equation  (3a)  in  (I5a)  we  obtain 

x* 
(16)   0  <  c^  -c^  =  I  k(x*,^)[w-2(e;cp)  -  w-2(4;c^)]  d^  . 

0 

Since  k  ^  0  it  follows  from  (15b)  that  the  integrand  in 

(16)  is  non-positive.   Hence  (I6)  yields  a  contradiction 

and  the  theorem  is  proved. 
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THEOREM  4.      lim     x„(c)    =   oo . 


OD 


0^ 


Proof:   From  (5a)  we  have,  w(x  (c),c)  =  S 
Xq  o 

=  c  -  J  k(xQ,^)w"^(^;c)  d?.   Since  w(^;c)  >  S  in  [0,x^], 
we  have  c-S  Jl  ■^2-  J   ^^(^0^^^  ^^'      Therefore  when  c  ->  oo 
the  integral  ->  oo  .   Since  the  integrand  is  bounded  this 
implies  that  x„  — >  oo  . 

THEOREM  5.   Xq(c)  is  continuous  for  c  ^  S. 

Proof:   From  Theorem  2  we  have  that  w  and  w'  are 
continuous  f\inctions  of  c.   Since  w'(x-.(c);c)  <  0,  applica- 
tion of  the  implicit  function  theorem  to  w(x„'c)  =  S  proves 
the  theorem. 
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4.   The  BoiJndary  Value  Problems. 

The  basic  theorem  for  Problem  S  follows  easily  from 
the  properties  of  the  initial  value  problem  established  in 
the  previous  section. 

THEOREM  S.   Problem  S  has  a  unique  solution  s(x). 
It  is  a  positive  and  monotone  decreasing  function. 

Proof.   Since  x^(c)  is  monotone  and  varies 
between  X  <  1  and  infinity,   there  is  a  \inique  value  Cq 
of  c  such  that  x^(c)  =  1.   The  positive  and  monotone 
decreasing  function  s(x)  =   w(xjCq)  is  the  imique  solution 
of  (2.7)  that  satisfies  the  boundary  condition  (2.6a)  and 
hence  is  the  required  solution  of  the  Problem  S. 

When  S  =  0  in  (2.6a)  Q/s   and  hence  s'  may  be  singular 
at  X  =  1  and  the  above  analysis  must  be  modified.   We  shall 
call  Problem  S  with  S  =  0,  Problem  0.   We  consider  a 
sequence  of  boundary  value  problems, 

(1)       Ls^  =  0  ,   s^(0)  =  0  ,   s^(l)  =  S^  >  0  ,   n  =  l,2,.., 


w 


here  {S  }  is  a  monotone  strictly  decreasing  sequence  of 


n 


positive  numbers  that  converges  to  zero.   It  follows  from 
Theorem  S  that  for  each  n,  (l)  has  a  unique  solution  s^(x) 


4f 

Here  X  is  the  value  defined  in  Section  3  immediately 
preceding  Theorem  3. 
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which  also  solves  the  integral  equation 

X 

(2)  s^(x)  =  c^  -  J  k(x,^)  s-2  d^  . 

0 

Here  c   is  the  uniquely  determined  initial  value  c  =  s  ( O) . 
n  ^       J  n   n^  ' 

In  the  Appendix  we  show  that  the  sequence  {s  (x)}  converges 
uniformly  on  [0,1]  to  the  function  s(x)  with  s(l)  =  0. 
The  limit  function  is  independent  of  the  choice  of  the 
sequence  CS  3.   Therefore  taking  limits  in  (2)  we  observe 
that  s(x)  satisfies  the  integral  equation 

X 


(5)  s  =  c  -  J  k(x,e)  s"^  d?  , 

0 


where  c  =  lim  c   and 
n-xx>   ^ 


X 


(4)  s' .  =   -  x^  J  Q  s  ^  d|  . 

0 

Since  s  satisfies  (5):?  it  is  a  solution  of  the  following 

boTondary  value  problem  which  we  call  Problem  B.   It 

consists  of  determining  functions  with  continuous  first 

derivatives  on  (0,1)  which  satisfy  (2.1)  on  (0,1)  and 

(2.5)  and  (2.6a)  with  S  =  0.   As  observed  by  Srubshchik  [6] 

a  solution  of  Problem  B  is  not  necessarily  a  solution  of 


If  p(r)  is  a  homogeneous  function  of  degree  m,  then  the 
existence  theorem  for  this  boundary  value  problem  can  be 
easily  established  by  considering  the  change  of  variables 
y=  c~^x,  G(y)  =  s(x)/c  where  a  =   3/(2m+5). 
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Prottlem  0  since  s'(l)  may  be  unbounded.    We  now  give  a 
necessary  and  sufficient  condition  for  Problem  0  to  have 
a  unique  solution. 

THEOREM  0.   Problem  0  has  a  unique  solution  if  and 
only  if  Q(l)  =  0. 

Proof:   Let  s(x)  be  the  unique  solution  of  Problem  0. 
Then  s(x)  is  also  a  solution  of  Problem  B  and  s'(x)  is 
continuous  on  [0,1].   Hence  from  (4),  -s'(l)  =  I 

=  I  Q,  s"^  d^   exists.   Assume  that  Q(l)  ^   0.   Since 

p      — P         ?       —1  ' 
I  exists,  then  so  does    s'  s   d^  =  -     (s   )   d^. 

^0  0 

However,  this  integral  is  unbounded  and  we  have  a 

contradiction.   Thus  Q(l)  =  0  is  a  necessary  condition 

for  Problem  0  to  have  a  unique  solution.   To  prove  the 

sufficiency  we  integrate  (4),  for  any  solution  s  of 

Problem  B,  from  x  to  1  and  use  the  boundary  condition 

s(l)  =  0  and  the  inequality  s(x)  <  c  to  obtain 

1     Ti 

(5)  s(x)  >  c"^  f(x)  =  c"^  I  '^'^  I  ^   ^^  ^T  • 

X     0 

The  integral  (5)  exists  since   lim  -^         \      Q,  d^  =  0. 


r\->Q)  0 

Thus  by  inserting  (5)  in  s'(l),  we  see  that  it  is  bounded 


We  recall,  see  Section  2,  that  Problem  0  requires  that  s' 
is  continuous  on  the  closed  unit  interval. 
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if  f  Q  f~  di   exists.   We  observe  from  (2.2)  that  Q,(l)  =  0 
implies  that  Q,'  (l)  =  0.   The  function  f  (x)  is  monotone 
decreasing,  f(l)  =  0  and  f'(l)  /  0.   Hence  Qf^  is 
bounded  at  x  =  1  and  Problem  0  has  a  solution.   The 
uniqueness  is  established  in  the  usual  way  by  assuming 
that  Problem  0  has  two  solutions  s,  and  Sp.   They  both 
satisfy  (3).   Since  s-|^(l)  =  s^Cl)  =  0,  we  have  from  (3) 

1 
(6)       s^(0)  -  s^CO)  =  J  k(l,^)  [s-2  -  s-2]  d^  . 

0 

From  Theorem  1  we  may  choose  s-,(o)  >  3^(0)  and  by  the 
proof  of  Theorem  ^,    s  (x)  >  s  (x)  for  0  <^  x  <  1.   Thus 

(6)  yields  a  contradiction  and  s-,  (x)  =   s  (x).   The 
condition  Q(l)  =  0  implies  that  the  pressure  has  a  zero 
resultant,  i.e.  p(r)  is  self-equilibrating. 

The  corresponding  result  for  Problem  H  is 
established  by  showing  that  for  the  function  g  defined  by 

(7)  g(x;c)  =   x^  w'(x;c)  +  aw(x;c)  , 


where  w  is  the  solution  of  the  initial  value  problem  (2.7), 
there  is  a  unique  value,  c„,  of  c  such  that 


(8)  g(i;cQ)  =  0. 
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Then  s(x)  =   w(xjc^)  is  the  lonique  solution  of  Problem  H. 
We  shall  now  obtain  the  necessary  properties  of  g. 

From  the  definition  (7)^  g  and  g'  are  obviously 
continuous  functions  in  [O^x^]  for  any  S  >  0.   Substitut- 
ing the  integral  relations  {j>.3)    and  {3-^)    in  (?)  we  obtain 

X 

(9)  g  =  ac  -  r  K(x,q  w"^(C)  dC  .  K{x,^)=  aK{x,^)+Qi 

d 

0 
Differentiation  of  (9)  gives 

X 

(10)  g'   =  -  QAj^  -  ax"^  J  (QA^)  d^   <   0  . 

0 

The  strict  inequality  holds  in  (lO)  if  Q  ^  0  in  (0,x). 
Thus  for  each  c,  g  is  a  monotone  decreasing  function  of  x. 
Therefore  if  c  <  0  it  cannot  have  a  positive  zero  and  we 
need  only  consider  c  >^  0. 

From  i'^.^)    and  (3.1)  we  observe  that  w'  <  0  and 
for  sufficiently  large  x,    w'  assiunes  arbitrarily  large 
negative  values.   Therefore  by  making  S  sufficiently  smallj 
we  conclude  from  (7)  that  g  has  a  zero  at  x  =  y(c)  in  the 
interval  0  <  y  <  x^.   The  strict  inequality  must  hold  in 
(10)  in  some  neighborhood  of  y.   Therefore  x  =  y  is  the 
unique  zero  of  g  for  fixed  c.   The  function  y(G)  has  the 
properties  summarized  in  Theorem  6. 
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THEOREM  6.   y(c)  Is  a  contlnuouSj  strictly 

monotone  Increasing  function  for  c  >  0  and  11m  y(c)  =  Y  <  1, 

c-^0 
11m  y(c)  =  oo  . 
c->oo 

The  proof  of  this  theorem  will  not  be  presented 
since  it  employs  arguments  similar  to  those  used  In 
establishing  Theorems  3,  4  and  5. 

Thus  from  Theorem  6  we  conclude  that  there  is  a 
unique  value  c„  of  c  for  which  y(c^)  =  1  and  hence 
g(l;c„)  =  0.   Therefore  we  have  established 

THEOREM  H.   Problem  H  has  a  unique  solution  s(x). 
It  is  a  positive  and  monotone  decreasing  function. 
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5.   Numerical  Solution. 

The  method  of  proof  used  in  the  previous  sections 
is  constructive.   We  shall  now  apply  it  to  obtain  a 
numerical  solution  of  Problem  H  with  p  =  constant.   From 
(2.2), 

Q(x)  =  Px^  ,      P  =  (l/8)(pR/Et)^  . 

We  define  a  new  dependent  variable  T  by 

T(x)  ^  s(x)V^  . 
Then  Problem  H  is  given  by 

(1)  (x\')'  =  -  x^/r^  ,     t'(o)  =  t'(i)  +  I^T(1)  =  0 

From  Theorem  H  in  the  previous  section,  it  follows  that 
there  exists  a  unique  value  of  c  =  T(0)  such  that 

(2)  f(c)  =  t'(i;c)  +  IJ.  T(l,c)  =  0  . 

We  determine  the  value  of  c  numerically  by  solving  (2) 
by  Newton's  method.   To  do  this  we  must  consider  the 
variational  equations  with  respect  to  c  of  the  initial 
value  problem.   We  have  solved  the  initial  value  problem 
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and  the  variational  problem  niimerically  by  the  Runge-Kutta 
method.    For  details  of  the  niomerical  application  of  the 

shooting  method  see  [11]-   The  values  of  T(x)  and 

1  A 
T  (x)  =  s  (x)/P  ' -^   obtained  from  this  computation  are 

graphed  in  Figure  1.   The  numerical  value  of  c  that  we 

obtained  is  0.8669531. 


In  the  calculations  we  used  the  value  v  =  0.52.   They  were 
performed  on  the  CDC  6600  Computer  of  the  AEC  Computing  and 
Applied  Mathematics  Center  of  the  Courant  Institute  of 
Mathematical  Sciences.   The  computation  was  done  by 
Dr.  L.  Bauer. 
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Appendix. 


THEOREM.   The  sequence  of  functions  {s  (x)3  defined 


n 


t>y  (^.1)  is  uniformly  convergent  on  [0,1]. 

Proof:   To  prove  the  theorem  we  must  first  show  that 

^^)  ^n+l^^)  ^  ^n^'')   ^°^   O^x^l.   n=l,2,... 


First  we  consider  x  =  0  where  s^(0)  =  c^  and  we  assume  thatt 
^n+1  "*  ^n-   ^^^  *^^^^  i^  a  point  x*  in  (0,1)  at  which 


^"^  ^n+1  ^   ^n     foi-  0  1  X  <  X*  .   Inserting  (4.2)  in  {2), 
we  obtain 

X 

^^)  =n+l-^n  =  Jk(x,?)  (s-2^  -  s-2)  d?  . 

0 

The  right  hand  side  of  (3)  is  negative.   This  provides  a 

contradiction  and  hence  c^  ,  <  c^ .   Since  c   >  0,  the 

n+±  n  n 

sequence  [c^]  converges.  By  a  similar  argument  we  can 
show  that  s^_^^  ^  s^  for  any  x  in  (0,1).  Therefore  (1) 
is  established.   We  now  consider  for  n  >  m  the  quantity 


I  We  need  not  consider  c   =  c^  ,  because  of  Theorem  1  and 

n+1   ^n 
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(s  -s  )'.   Using  (4. 2)  J  it  is  given  by 

X 

^   0 

Since  (l)  is  true  for  all  n,  the  right  hand  side  of  (4) 

is  positive  and  s  -  s   is  monotone  increasing.   Therefore 
■^  m   n 

we  have 


(5)       s^(x)  -  s^(x)  <  S^-  S^   for  0  <  X  <  1  . 


Since  the  sequence  [S  3  converges  to  zero  the  sequence 

{s  (x)3  is  uniformly  convergent  and  the  theorem  is  proved, 
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